Electron-phonon vertex in the two-dimensional one-band Hubbard model 



Z. B. Huang, W. Hanke, and E. Arrigoni 
Institut fur Theoretische Physik, Universitdt Wiirzburg, am Hubland, 97074 Wurzburg, Germany 

D. J. Scalapino 

Department of Physics, University of California, Santa Barbara, California 93106-9530 USA 

(Dated: February 2, 2008) 

Using quantum Monte Carlo techniques, we study the effects of electronic correlations on the 
effective electron-phonon (el-ph) coupling in a two-dimensional one-band Hubbard model. We con- 
sider a momentum-independent bare ionic el-ph coupling. In the weak- and intermediate-correlation 
regimes, we find that the on-site Coulomb interaction U acts to effectively suppress the ionic el- 
ph coupling at all electron- and phonon- momenta. In this regime, our numerical simulations are 
in good agreement with the results of perturbation theory to order U 2 . However, entering the 
strong-correlation regime, we find that the forward scattering process stops decreasing and begins 
to substantially increase as a function of U, leading to an effective el-ph coupling which is peaked 
in the forward direction. Whereas at weak and intermediate Coulomb interactions, screening is 
the dominant correlation effect suppressing the el-ph coupling, at larger U values irreducible vertex 
corrections become more important and give rise to this increase. These vertex corrections depend 
crucially on the renormalized electronic structure of the strongly correlated system. 



The role of the el-ph interaction in the physics of the 
high T c cuprate superconductors remains unclear. On the 
one hand, the linear T dependence of the resistivity up 
to high temperatures and the small value of the isotope 
coefficient of the optimally doped materials suggest that 
the el-ph interaction plays a secondary role 1 . The fact 
that the undoped cuprates are Mott antiferromagnetic 
insulators supports the notion that strong Coulomb in- 
teractions are dominant and that the essential physics is 
contained in the Hubbard and t-J models^. On the other 
hand, however, a variety of experiments also display pro- 
nounced phonon and electron-lattice effects in these ma- 
terials: superconductivity-induced phonon renormaliza- 
tion^, large isotope coefficients away from optimal dop- 
ingi, tunneling phonon structures^, etc., give evidence of 
significant el-ph coupling. Recently, photoemission data 
indicated a sudden change in the electron dispersion near 
a characteristic energy scaled, which is possibly caused by 
coupling of electronic quasiparticles to phonon modes. 

To elucidate the effects of strong electronic correla- 
tions on the el-ph interaction, several authors have calcu- 
lated the el-ph vertex function in the one- and three-band 
Hubbard models based on (1/iV) expansion within slave- 
bosonL& and X operator— formalisms. One finding^ is 
that for the ionic, i.e. onsite, el-ph coupling in the under- 
doped regime, the backward scattering with large phonon 
momentum transfer is suppressed much more than the 
forward scattering with small phonon momentum trans- 
fer. Based on this finding, it was argued that a forward 
peaking of the renormalized el-ph vertex could account 
for the absence of phonon features in the transport data. 
In addition, an el-ph interaction which is peaked at small 
phonon momentum transfer contributes to an attractive 
interaction in the d :c 2_ !y 2-pairing channel2*ifi. One should 
notice that these previous calculations were limited by 
the approximate nature of 1/N and slave-boson treat- 
ments, moreover, they were carried out for U — > oo. 



Here, we would like to use a more accurate numeri- 
cal method to gain further insight into the way in which 
strong electron correlations dress the el-ph coupling for 
a range of values of U from weak to strong correlation. 
This analysis is important since it turns out that the 
U dependence of the el-ph coupling shows a dramatic 
change as a function of U (see below). Specifically, we 
apply the determinantal Monte Carlo 11 algorithm to find 
the single-particle response to external phonon fields in 
the one-band Hubbard model. In particular, we will cal- 
culate an effective el-ph coupling 12 g(p, q) (effective el-ph 
vertex) for scattering quasiparticles near the Fermi sur- 
face (which includes screening, vertex corrections, and 
the quasiparticle renormalization) produced by the Hub- 
bard U. 

Our principle findings are: (1) Initially (up to the value 
of U w 6i), as the Hubbard- Coulomb interaction U in- 
creases, the el-ph coupling is suppressed by electronic cor- 
relations for all phonon and electron momenta, and in 
particular for the backward scattering around q = (jr, tt). 
The suppression is due to the conventional screening 
term. (2) The behavior changes even qualitatively in the 
strong-correlation regime (U > 6t). Here, the effective 
el-ph coupling at small phonon momentum transfer in- 
creases with increasing U, while the one at large phonon 
momentum transfer appears to saturate (see Fig.|3J. The 
increase of the el-ph coupling in the forward direction is 
due to irreducible vertex corrections, which become the 
dominant correlation effects at strong Coulomb interac- 
tions. These vertex corrections are intimately connected 
with the renormalized electronic structure. It is argued 
that the picture of a "spin-bag" quasiparticle can explain 
the qualitatively different el-ph couplings for small and 
large phonon momenta. Furthermore, we would like to 
stress that our numerical results for the charge compress- 
ibility, which decreases monotonically with increasing U, 
rule out the explanation of the increase of the el-ph ver- 
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tex at small q as a function of U in terms of a close-by 
phase separation or charge instabilitjii^. 

Our starting point is the one-band Hubbard model, 
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The operators c\ a and c ia as usual create and destroy 
an electron with spin a at site i, respectively and the 
sum (ij) is over nearest-neighbor lattice sites. Here, U 
is the onsite Coulomb interaction and we will choose the 
nearest-neighbor hopping t as the unit of energy. 

In our simulations, we have used the linear-response 
technique in order to extract the el-ph vertex function. 
In this method, one formally adds to Eq. the interac- 
tion with a momentum- and (imaginary) time-dependent 
lattice distortion (phonon) field Uqe~ lqaT in the formi2ii4 
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where g® q is the bare el-ph coupling. One then considers 
the "anomalous" single-particle propagator in the pres- 
ence of this perturbation defined asi£ 
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where () H+H ct ^ ph is Green's function evaluated with the 
Hamiltonian H+H e i^ p h- Diagrammatically Ga{p, <?) has 
the structure shown in Fig. ^ so that the el-ph vertex 
function T(p, q) can be expressed quite generally in terms 
of Ga and of the single-particle Green's function G(p) in 
the form 



T(p, q) = lim 
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It is, thus, sufficient to calculate the leading linear re- 
sponse of Ga to H el-ph j which is given by 

G A (P, q) = U n f dre^+^r j dT ' e ^or' J2 g q x 
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where + is a positive infinitesimal. The two-particle 
Green's function in Eq. (0 is evaluated with respect to 
the pure Hubbard Hamiltonian (Eq. QJ). Since we have 
only considered the linear-response contribution from the 
phonon field, the el-ph vertex T contains full contribu- 
tions from Coulomb interactions onlyi£, 

Close to the Fermi energy, the single-particle Green's 
function can be written as 



G( P ) 
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where Z(p) is the wave-function renormalization and E p 
the quasiparticle excitation energy^. Then for electron 




G A (p,q)= 



FIG. 1: Diagrammatic representation of Ga(p,q) within lin- 
ear response to it q . The thick solid lines represent dressed 
single-particle Green's functions of the Hubbard model. The 
wavy line denotes the external perturbation in Eq. 
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FIG. 2: Low-order Feynman diagrams for the irreducible el-ph 
vertex A(p, q) (top) and low-order polarization graphs (lower) 
that enter the full vertex F. The thin solid lines are the non- 
interacting Green's functions and the dashed lines represent 
the Hubbard interaction U. The thin (thick) wavy line stand 
for the bare (screened phonon) fields. 



scattering processes which involve states near the Fermi 
surface, the effective el-ph coupling reads 
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In the following, we will focus on the case of an ionic 
el-ph coupling, in which the bare coupling (/p q is a con- 
stant g°. Since we are considering linear terms in g° only, 
we can set g° equal to 1. This corresponds to the sim- 
ple Holstcin form of the el-ph interaction, which is an 
important limiting case. Moreover, having the bare in- 
teraction independent of p and q makes it easier to see 
modifications, which arise from the strong U correlation 
effects. 



The low order U and U 2 vertex contributions to T are 
displayed in Fig. [3 The diagrams shown at the bottom 
of Fig. |21 are the leading terms of the RPA approxima- 
tion (1 — i U n ((7)) _1 to the polarization correction, with 
Ho{q) the contribution from the single bubble. The ex- 
act Monte Carlo result for the polarization correction is 
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With this in mind, T can be written in terms of the 
screening factor and an irreducible vertex A(p, g), which 
is the sum of graphs that can not be separated into two 
pieces by cutting a single dashed Coulomb interaction 
line U (see Fig. |5J), i-e. 



T(p,q) = {l + \im{q))k(p,q) 



(9) 



Thus, the strong-correlation effects associated with the 
Hubbard-Coulomb interaction U lead to an effective el- 
ph coupling which can be expressed in the canonical form 



{l + ^UU(q))A(p,q) 
(Z(p)Z(p + q))h 
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and one sees that it consists of a product of an irreducible 
vertex A(p,q), a screening factor (1 + |?7II(g)), and a 

quasiparticle renormalization (Z(p)Z(p + <z))~ 5 factor. 

Our numerical Monte Carlo simulations were per- 
formed on an 8 x 8 lattice at an inverse temperature (3 = 2 
and a filling < n >— 0.88. We have set the frequencies 
to their minimum values, i.e., pa — ttT for fermions and 
go = for bosons. We have checked some special cases 
for which one can reach lower temperatures, namely a 
2D system at weak correlation and/or with large doping 
(< n >= 0.65). For these systems, we have found that 
the real part of the vertex function r((7rT, p), (0, q)) de- 
pends only weakly on temperature, and the imaginary 
part always vanishes as T — > 0. In the following we will, 
therefore, focus on the real part of the vertex function 
at po — 7rT. Comparison with exact diagonalization on 
a 4-site ring demonstrates that the difference in T(p, q) 
between the two results is less than two percent up to 
U = 8. 

We are interested in el-ph scattering processes in which 
the incoming and the outgoing electron momenta p and 
p + q are close to the Fermi surface. For an 8 x 8 lat- 
tice doped near half-filling, the q and U dependence of 
g(p, q) for the scattering processes on the half-filled di- 
amond Fermi surface is studied. In particular, we will 
examine initial states corresponding to p = (— 7r, 0) and 
p = (— 7r/2, 7r/2). Other choices of p and p + q close to 
the half-filled diamond Fermi surface give qualitatively 
similar results to those reported here. 



Monte Carlo results for g(p, q) and for the polariza- 
tion factor (1 + ^£7n) are shown in Fig. [21 The left 
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FIG. 3: Real part of the effective el-ph coupling g(p,q) and 
the polarization factor 1 + \UYl versus q for (a) U < 6 and 
(b) U > 6. Here q = Tr(h, h) with h the tick label of the x 
axis. The incoming electron carries momentum p = (— n, 0) 
and the value of U is indicated by the shape of the symbol. 



hand side, Fig. [2Jl, shows the behavior in the weak- and 
intermediate-correlation regimes. The right hand side of 
the figure, Fig.[2p, shows similar results when the system 
enters the strong-correlation regime. One can clearly see 
that when the Hubbard U is smaller than U ~ 6, g(p, q) 
decreases as a function of U from its bare value g° = 1, 
for all momentum transfers. Then, as the interaction U 
increases to the strong-correlation case (U ~ W = 8t), 
the effective el-ph coupling begins to increase. This be- 
havior is particularly evident at smaller values of mo- 
mentum transfer. Our finding at large phonon momen- 
tum is similar to that of Deppeler et aZ.'s workii, which 
shows that the local el-ph interaction is suppressed by 
both electronic correlations and dynamic phonon vertex 
corrections. In the strong-correlation regime, the over- 
all q dependence of the el-ph coupling agrees reason- 
ably well with the results of the 1/N expansion^, which 
are obtained for the U — > oo limit. However, in our 
case the interesting behavior is that the effective el-ph 
coupling as a function of U is nonmonotonic, first de- 
creasing and then, at physically interesting values of U, 
increasing. This finding deviates from the prediction 
of a Fermi-liquid analysis 8 . According to this analysis, 
lim q _ ) ,o,go=0 9{Pi ?) u F » with Fq the zero-harmonic 
symmetric Landau amplitude so that the effective el-ph 
coupling decreases monotonically with increasing U since 



Fq becomes larger with U (except when approaching a 
charge instability). 

From Fig. [21 one can see that the polarization factor 
acts quite generally to suppress the el-ph coupling. At 
large momentum transfer, this quantity saturates as U 



4 




2 4 6 8 10 

u 

FIG. 4: Real part of g(p, q) as a function of U for (a) p = 
(—it, 0) and (b) p = (— 7r/2, 7r/2). The value of q is indicated 
by the shape of the symbol. The solid circles are Monte Carlo 
results and the open symbols show the perturbation theory 
contributions shown in Fig. |5| 

increases, while at small momentum transfer it continues 
to decrease. Comparison between g(p, q) and the polar- 
ization factor indicates that at weak correlation screening 
is the dominant correlation effect suppressing the el-ph 
coupling. On the other hand, with an increasing Hub- 
bard U, the vertex corrections A become more impor- 
tant, making the effective el-ph coupling peaked in the 
forward scattering direction. 



In order to see the U dependence more clearly, in Fig.|3| 
Quantum Monte Carlo calculations are compared with 
perturbation theory for different values of U. Here, the 
solid symbols are Monte Carlo results and the open sym- 
bols show the results obtained by evaluating T(p, q) per- 
turbatively with the diagrams of Fig. [21 In the pertur- 
bative calculations, g(p, q) is calculated by using wave- 
function renormalizations Z(p) and Z(p + q) extracted 
from Monte Carlo data. As one can see, in the weak- 
correlation regime, the perturbative calculations are in 
good agreement with Monte Carlo simulations. However, 
when the Hubbard U exceeds U ~ 4 (~ perturbation 
theory appears to break down. 

QMC calculations of the doping, temperature and U 
dependence of the vertex enhancement and its q depen- 
dence, which will be presented in more detail in a longer 
paper, give a physical picture for both the unexpected 
increase of the vertex as a function of larger [/-values 
for small phonon momenta and also the suppression of 
the vertex for large phonon momenta: It is well estab- 



lished, in particular in terms of QMC work on the single- 
particle spectral function of the Hubbard modeli&, that 
the single-particle excitations as a function of increas- 
ing [/-values undergo (at lower enough temperatures, 
around f3 = 2 — 3) a crucial physical transition into a 
strong-correlation regime around U ~ 6t: the electroni- 
cally filled valence band of width W, which is essentially 
given by the bare band width W = 8t, splits into two 
"bands". The physical picture behind this splitting is 
the formation of a "spin-bag" quasiparticle, i.e. of the 
bare particle (hole) dressed with a spatially (typically a 
few lattice constants) extended spin cloud, which is due 
to the frustration of the local antiferromagnetic order. 
The spin bag moves coherently and "slowly" with an en- 
ergy scale J = within the new, strongly renormal- 
ized quasiparticle band of width ~ J. This coherent mo- 
tion couples effectively (with small energy denominators) 
to longer wavelength lattice displacements, whose wave- 
length is typically longer than the spin-bag "diameter" . 
On the other hand, there is at larger [/-values also an in- 
coherent lower Hubbard band whose higher energy scale 
corresponds to the "rattling around" of the bare parti- 
cle within the spin bag«i&. Only large momenta phonons, 
i.e. with wavelength smaller than the "extension" of the 
spin bag, can couple to these incoherent electronic de- 
grees of freedom. Their coupling is weak because of the 
combined effects of the incoherent motion and the large 
(~ scale W) energy denominators. 

In summary, based on QMC simulations, we have stud- 
ied the el-ph vertex function in the two-dimensional Hub- 
bard model. We find that in the weak-correlation regime, 
the effects of the Hubbard interaction U are to suppress 
the ionic el-ph coupling at all phonon momenta, with 
backward scattering processes being more strongly sup- 
pressed than forward ones. On the other hand, in the 
strong-correlation regime, the vertex at smaller phonon 
momentum transfer anomalously increases as a function 
of U. We also find that screening is the dominant con- 
tribution to the vertex corrections at weak correlation, 
while at strong correlation the irreducible vertex correc- 
tions are crucial. 
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